CLAIMS 

1.-52. (CANCELLED) 

53. (CURRENTLY AMENDED) A method of estimating arterial delay and arterial 
dispersion (t, a, a) values for outputting blood perfusion indices for a region of 
interest (ROI) by from operating a computer program on intensity data in input to 
a computer comprising: 

a. using a computer to apply applying a first gamma-variate function 
(GVF) to an arterial input function (AIFa) using a computer t o provide an 
estimated first model of a vascular transport function ha(t), wherein for t < 

ti, ha(t) = 0 and for t > ti, ha(t) = —(t-t^Y^e'^^'^^^^^^ , wherein m 

e stimat e d ti is the transit time of a contrast agent from a measured initial 
said AIFa to a region of interest (ROD and 0 ] _ is an estimating an initial 
delay value of said contrast agent, wherein said Oi=(t } )(^ \ )/(l'^ \ ). wherein 
said is a known relative dispersion value having a range from 0 to 1: 

h: e stimating an initial valu e 04 . of said contrast ag e nt, wh e r e in said 

e4.=(t^)(^4 ^)/(l -#4 -), wherein said § 4- is a known relative disp e rsion value 
having a range from 0 to 1 ; 

b. using a computer to convolve convolving A IFaftl with said hft(t7^=0) 
haffl with 0^=0 using a computer for obtaining an arterial input function 
AIFt(t) = AIFa(t) ® hft(treM.=0) hJt) with a^=0 at said ROI; 

c. using a computer to estimate estimating a blood flow rate Ft and a tissue 
impulse residue function R€(t) using a computer b v deconvolving a 
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concentration curve C(t) = (Ft/kH)AIFt(t) ® Re(t), wherein kn is a 
hermocri t hermatocrit correction constant having a known value; mi 

d. using a computer to optimize said mean transit time and dispersion (t7. 

values using a least squares method: and 

e. using a computer to output said outputting estimated and optimized 
tissue mean transit time and dispersion (t2, 02, 02) values from an 
estimated transport function he(t) for input to a simulated transport 
function hs(t), wherein a simulated tissue impulse residue function Rs(t) is 
determined, wherein a simulated concentration curve Cs(t) is fitted to said 
measured C(t) and quantitative said blood perfusion indices are calculated^ 
wherein each said step is performed bv a suitablv programmed computer . 

54. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said 
intensity data is generated by administering a contrast agent to a body lumen 
of a body during a dynamic imaging scan, wherein said body Ixmien comprises 
an artery or vein, wherein an image response from said contrast agent is 
recorded to computer data storage in a computer. 

55. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said C(t) is a 
temporal concentration of said contrast agent obtained from said intensity data, 
wherein said intensity data comprises contrast images sequentially acquired 
from a region in a body, whereby said contrast agent concentration is plotted 
versus time. 
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56. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said AIFa is 
based on a measured early arrival contrast agent peak intensity from a feeding 
blood vessel to said ROI. 

57. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said AIFa is 
scaled upward according to a venous input function (VIF), wherein said VIF is 
based on a measured late arrival contrast agent peak intensity from a large vein 
draining from said ROI. 

58. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said 
estimated transit time ti is the transit time of said contrast agent from a 
measured initial said AIFa of said contrast agent C(t) in a body lumen to said 
ROI, wherein said ti is estimated from plots of said AIFa versus time and said 
C(t) versus time. 

59. (CURRENTLY AMENDED) The method of claim 53, wherein said ha(t) is 
calculated using said estimated transit time ti and said estimated dispersion 
value Gu wherein h^ftr-ea=Q 4 hJt) with a2^=0 is plotted versus time. 

60. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said 
estimated transport function he(t) is calculated using the relation he(t) = - 
dRe(t)/dt. 
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61. (CURRENTLY AMENDED) The method of claim 53, wherein said tissue 
mean transit time and dispersion (t2, aa, 02) values are estimated from said 
estimated transport function he(t), wherein said t2, said 02 and said a2 are input 
to a simulated transport function hs(t), wherein said hs(t) is said a second 
gamma- variate function. 

62. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said 
simulated tissue impulse resistive function Rs(t) is determined using the 



63. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said 
simulated concentration curve Cs(t) is determined using the relation Cs(t) = 



64. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said Ft, said 
ti, said Oi, said a2, said t2, said a2, and said 02 are optimized by a least squares 
method to fit said Cs(t) to said C(t). 



65. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said 
perfusion indices have the relations: 
f blood flow (BF) = Ft; 




0 



(F,/kH)AIF,(t) ® Re(t) = (FAh) / AIF,(t) R,(t-T)dT. 
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g. Mean Transit Time (MTT) = t2 + 02(1 +a2); 

h. Blood Volume (BV) = BF * MTT; 

i. Arterial Delay (DT) = ti + ai(l+ai); 

j . Arterial Dispersion time (ADT) = 0\ Vl + oci ; 

5 k. Tissue Dispersion Time (TDT) = 02 Vl + aa ; 

L Relative Arterial Dispersion (RAD) = ADT/DT; and 
m. Relative Tissue Dispersion (RTD) - TDT/MTT. 

66. (PREVIOUSLY PPIESENTED) The method of claim 53, wherein said AIFt(t) 
is measureable in a small lumen showing a delay relative to said AIFa(t), 
wherein optimized values for said 0\ and said ti are determined by fitting said 
simulated AIFt(t) to said measured AIFt(t), wherein said relative dispersion Pi 
is determined and applied to all other said intensity data of said ROI using said 
Pi, wherein a more robust fitting process is provided by a reduced number of 
parameters for optimization. 

67. (PREVIOUSLY PRESENTED) The method of claim 66, wherein when 
said relative dispersion Pi is determined, said vascular transport fiinction 
ha(t) is described by a single variable said ti with a constant said Pi, 
20 wherein a two-step method is used to determine said delay and said 

dispersion values comprising: 

a. deriving an initial tissue impulse residue fiinction Ro(t) by 
deconvolving C(t) = (Fo/kH)AIFa(t) ® Ro(t) using a model-free singular 
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value decomposition (SVD) method, wherein said time delay ti is 
determined by a maximum position of said Ro(t) at Ro max(t=ti); and 
b. determine said AIFt(t) at an input of said ROI using said ha(t) with said 
ti and said Pi held constant, wherein said 0\ is determined. 

68. (PREVIOUSLY PRESENTED) The method of claim 67, wherein a 
value of tissue blood flow Ft and a corrected impulse residue function 
Re(t) are obtained by deconvolving C(t) = (Ft/kH)AIFt(t) ® Re(t) using 
said SVD method, wherein said perfusion indices are determined from 

a curve of said R^Ct), wherein MTT= Rc(T)dT, BF=Ft, and 




BV=BF*MTT. 

69. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said contrast 
agent is in a tissue ROI having a tissue mean transit time r, wherein a tissue 
impulse residue function is approximated by the relation R(t >t) = Ee"^*"^^ and 
R(t^T) = 1, wherein E is an extraction fraction of said contrast agent in said 
tissue, wherein k is a constant clearance rate of said contrast agent diffusing 
from said tissue having a relation k = E*FtA^e, wherein Vgis the volume 
fraction of extravascular and extracellular space (EES) in said tissue. 

70. (PREVIOUSLY PRESENTED) The method of claim 69, wherein said 
tissue impulse residue function Rs(t) of said simulated concentration curve 
Cs(t) is replaced by an average impulse residue function that incorporates 
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said contrast agent leaked out of a blood vessel into said tissue and 
gradually clearing from said tissue, wherein said simulated concentration 
curve Cs(t) is fitted to said measured C(t) and quantitative said blood 
perfusion indices are calculated, wherein said E and said Ve are additional 
5 parameters optimized with other adjustable parameters using a least squares 

method. 
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